Tunable Casimir repulsion with three dimensional topological insulators 
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In this Letter, we show that switching between repulsive and attractive Casimir forces by means 
of external tunable parameters could be realized with two topological insulator plates. We find two 
regimes where a repulsive (attractive) force is found at small (large) distances between the plates, 
canceling out at a critical distance. For a frequency range where the effective electromagnetic action 
is valid, this distance appears at length scales corresponding to 1 — e(aj) ~ 



The full experimental accessibility to micrometer and 
sub-micrometer size physics and the possibility of devel- 
oping applications has turned the understanding of phe- 
nomena at these scales to be of fundamental importance. 
Within this scenario, the Casimir force [l| arises when 
two objects are placed near each other at distances of a 
few micrometers. In the general case of two dielectrics 
the situation is well described by the theory developed by 
Dzyaloshinskii et. al where the optical response of the 
material determines the magnitude and behaviour of the 
force. In the simplest case of a mirror symmetric situa- 
tion a theorem ensures that the force is always attractive 
SSl) resulting in a problem for nano- mechanical devices. 
To revert the sign, one must search non symmetric sit- 
uations, usually adding complexity to the problem. The 
first Casimir repulsion proposal, known as Dzyaloshinskii 
repulsion, was recently confirmed experimentally Q and 
it involves a third dielectric medium between the plates, 
excluding the possibility of frictionless devices and quan- 
tum levitation. In turn, vacuum mediated proposals in- 
clude magnetic versus non-magnetic situations [6| and 
the use of metamaterials 0-01 • In this Letter, we report 
a new method for obtaining a twofold tunable Casimir 
repulsion. By use of the optical properties of topologi- 
cal insulators (TI) it is feasible to achieve all situations 
between repulsion to attraction by using two control- 
lable parameters: the distance between dielectrics, and 
the sign of the topological magnetoelectric polarizability 
(TMEP) 9, where the latter allows to tune the optical 
properties of the mentioned materials. 
TIs are characterized by a bulk insulating behaviour with 
metallic boundary states protected by time reversal sym- 
metry 13, HI- The topological protection of edge states 
ensures their stability against non-magnetic perturba- 
tions. The 3D counterpart of this novel topological state 
was shown to exist in a Bi^^Sbi-a: alloy [l2l | and in the sto- 
ichiometric crystals Bi2Se3, Bi2Te3,TlBiSe2 and Sb2Te3 

The Casimir force is intimately related to the optical 
properties of the two dielectric bodies Q ■ For instance, 
consider the situation where two dielectric parallel semi- 
infinite bodies (labeled 1 and 2) are placed at a distance 
d from each other in vacuum. In this case, the Casimir 
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FIG. 1: Different configurations of the Casimir effect with 
identical TI covered with a thin magnetic layer. (a) The mag- 
netization is of the same sign on each surface resulting in a 
case where 9i = 62 giving Casimir attraction. In (b), (c) 
and (d) the magnetizations have opposite signs on the sur- 
face (^1 = —62) leading to attraction when d > repulsion 
when d < dm and to a quantum levitation configuration at 
dm where the net force is zero. 



energy density (CED) stored by the plates is given by 
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where A is the plate area, ks = ^fe^ -|- ^^/c^ is the wave 

vector perpendicular to the plates, A;|| is the vector paral- 
lel to the plates and ^ is the imaginary frequency defined 
a.s Lu — The matrices R1.2 are 2x2 reflection matri- 
ces of media 1 and 2 containning the Fresncl coefficients 
defined as: 



R 



i?s,s(i^,k||) i?5,p(i^,k| 
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where Ri^ describes the reflection amplitude of an in- 
cident wave with polarization i which is reflected with 
polarization j. The polarizations Rs (p) describe paral- 
lel (perpendicular) polarization with respect to the plane 
of incidence. The Casimir force per unit area on the 
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plates is obtained by differentiating expression ([Ij. A 
positive (negative) force, or equivalently a positive (neg- 
ative) slope of Ec{d), corresponds to attraction (repul- 
sion) of the plates. 

The electromagnetic response of a dielectric, which de- 
fines the reflection matrices, is governed by Maxwell's 
equations derived from the ordinary electromagnetic ac- 
tion 5*0 = / dx^dt {eE'^ - (l/^)^^), being E and B the 
electric and magnetic fields respectively. TI in three di- 
mensions are well described by adding a term of the form 
Se = (a/47r2) / dx^dt E ■ B, where a = 1/137 is the 
fine structure constant and 6 is the TMEP (axion field) 

11 , l3| . Because of time reversal symmetry, this term is 
a good description of the bulk of a trivial insulator (e.g. 
vacuum) when 6 = and of the bulk of a TI when 9 = tt. 
However, the axion coupling is only a good description 
of both the bulk and the boundary of a TI when a time 
reversal breaking perturbation is induced on the surface 
and the system becomes fully gapped. In this situation, 
9 can be shown to be quantized in odd integer values of 
TT such that 9 = {2n + l)7r where n G Z. The value of n 
is determined by the nature of the time reversal breaking 
perturbation, which could be controllable experimentally 
by covering the TI with a thin magnetic layer. In particu- 
lar, positive or negative values of 9 are related to different 
signs of the magnetization on the surface [l^ . As we will 
demonstrate in what follows, the Casimir force is very 
sensitive to the value of 9 and the tunability of its sign 
will allow us to describe a mechanism for switching be- 
tween repulsive and attractive forces. 

The electromagnetic response of a system in the pres- 
ence of a 6* term is still described by the ordinary Maxwell 
equations but the constituent relations which define the 
electric displacement D and the magnetic field H acquire 
an extra term proportional to 9 [19:] D = eE + a{9/TT)B 
and H = B / ji — (\{9/-k)E. We note that eq.(l) can be 
easily modified to take into account magnetoelectric cou- 
plings (this happens also in chiral metamaterialsL9]). The 
result is the same equation with the proper refiection ma- 
trices. It is then possible to derive by means of ordinary 
electromagnetic theory the reflection coeflicients of a TI- 
vacuum interface. For a TI characterized by a frequency 
dependent dielectric function e(a;) and a TMEP 9 the re- 
fiection coefficients will take a symmetric form [1] where 
the off diagonal coefficient can be expressed as: 



-Rs,p(jf,k||,6') = sgn(0)r^p(i^,k||,( 
where r^p is an even function of 9. When 9 - 
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0, Rs,p — 0, 

leading to the usual attractive Casimir force due to the 
non- mixing of polarizations [21|J. When 9 the reflec- 
tion coefficients mix polarizations and the sign of 9 plays 
a crucial role on the sign of the Casimir force (see below) . 
In what follows we will consider that the surface of the 
TIs of the Casimir system are covered by a thin mag- 
netic layer as shown in Fig.l. This effectively turns the 
TI into a full insulator (both in the bulk and on the sur- 



face) which can be safely described with the TMEP and 
a dielectric function, as shown in earlier works [III [3l ■ 
Hence, to numerically compute the CED by means of ([T]) 
a model for the dielectric function is necessary (we hence- 
forth assume /i = 1). Due to the low concentration of free 
carriers in insulators the most general phenomenological 
model to describe the optical response of a dielectric is 
a sum of oscillators to account for particular absorption 
resonances |2lJ. When only one oscillator is considered 
(see however [22|) the dielectric function evaluated can 
be written as: 



e{iO = 1 + 



-UJt 



(4) 



In this model, ujr is the resonant frequency of the os- 
cillator while We accounts for the oscillator strength. 
The damping parameter satisfies << uir play- 
ing therefore a secondary role. In what follows, we have 
rescaled all quantities in units of cor leaving the quan- 
tity e(0) = 1 -|- (we/wfl)^ as the only parameter of the 
model. A good candidate to be described by this model 
is the TI TlBiSe2[lB|- This material has experimentally 
[H (neglecting free carrier contributions and assuming 
high frequency transparency) e(0) ~ 4 and has a single 
resonant frequency near 56cm^^. Other TI could need 
more oscillators to be added in (|4]). 

We have computed the CED between two TI plates de- 
scribed by the TMEP 9i and 92 and the value of the 
dielectric constant at zero frequency e(0) by numerical 
evaluation of expression ^ . The results are summarized 
in figures [2] and [3] where the CED is plotted against the 
dimensionless distance d = duji/c. From Fig. [2] (a) it 
is clear that opposite signs of 6*1 _ 2 lead to the existence 
of a minimum (dm) where the net force is zero. The 
behaviour is attractive when both signs become equal 
suggesting that it is possible to tune the Casimir force 
by tuning the relative signs of 0, i.e switching the mag- 
netizations of the coverings. The existence of a mini- 
mum is analytically shown below in terms of the rela- 
tive importance of the off-diagonal terms (O against the 
diagonal terms and in terms of the relative sign of the 
TMEPs. At large distances the diagonal terms dominate 
and the usual Casimir attraction is recovered. At small 
distances, the off-diagonal terms dominate and their sign 
determines whether the CED approaches ±00 (i.e. repul- 
sive or attractive) leading to a minimum at intermediate 
distances if the signs of 9i^2 are opposite. 
To prove the existence of the minimum we consider the 
Fresnel equations for TI obtained earlier in |l| which 
lead to eq. ^ added to the following properties of the 
dielectric function: finite dielectric permittivity at zero 
frequency (e(0) < 00) and high frequency transparency, 
e(w) — > 1 when w — ?> cx). For analytical traceability we 
assume 9i = — ^2 = and that the dielectric function 
(HI describes the TI, although the derivation does not 
depend on the explicit form of the dielectric function as 
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FIG. 2: Casimir energy density (in units of £"0 ~ Ahc/ {lur/ c)'^) as a function of the dimensionless distance d for lOeJiOR — 
0.45. In (a) 6^2 = — tt is fixed. Whenever sgn(Si) = — sgn(^2) a minimum dm appears leading to a vanishing net force on the 
plates. Increasing Q\ within positive values suppresses the minimum shifting d„i towards lower values (if both signs are equal 
then only attractive behaviour occurs). Complete repulsion is achieved when one of the TMEP is much bigger than the other, 
(b) The optimal situation Q\ = — = d. Different values of 9 show that the minimum is enhanced when the difference between 
the two values is as small as possible {9 = n). 



long as it fulfils the mentioned conditions. 
In ([T]) we can rescale ^ and k|| to contain d which gives an 
overall factor 1 /cF and forces the reflection matrices to be 
evaluated at the rescaled frequency and momenta ^/d and 
k||/J. Hence, Ec{d — >• 0) — J- ±00 and Ec{d 00) ^> 
since the behaviour of e{i(,) ensures that the reflection 
matrices are not singular when evaluated at i^/d 
and i£,/d — >■ 00. 

The way the integral approaches these limits determines 
the sign of Ec{d). For instance, if the integrand is pos- 
itive at small distances and negative at large distances, 
necessarily a minimum exists at an intermediate distance 
dm- In what follows it will be shown that this is exactly 
what happens unless e(0) = 1, where both limits are pos- 
itive and hence long range repulsion is obtained. Under 
these conditions the diagonal terms in ([2|) are equal for 
both TI (which we label and rp) and the off-diagonal 
terms given by ([3]) have opposite overall signs, but equal 
absolute value given by the function r^p. Introducing 
these inside ([T]) the integrand reads: 

\og[l + e-^^i^\2r%~rl~rl)+ (5) 

where fcg is now evaluated at the rescaled frequency and 
momenta just as the reflection matrices. In the limit of 
small distances {d — >■ 0) and using the high frequency 
transparency of the dielectric function it can be shown 
that |rp| << \rsp\ since the first are of order and 
the second are of order a. Hence the integrand is positive 
and so Ec{d 0) — >■ -l-oo. 

Now we consider the limit of large distances {d 00). 
In this limit, the reflection coefficients take the form 
r, = (l- e(0) - a^) /D (a similar expression holds for 



rp) and r^p = 2|a|/£) where D = 1 + e(0) +a'^ + v/e(0)x 
, X is a frequency and momentum dependent function ir- 
relevant for the present discussion and a = aO/ir. In this 
long distance limit, depending on the values of e(0) dif- 
ferent behaviors emerge. Since e(0) > 1 we now consider 
the two extreme limits, one where e(0) = 1 and the other 
with e(0) >> 1. 

In the limit where e(0) >> 1, the condition \rs\, \rp\ » 
\rgp\ is always satisfied. When e(0) is strictly infinity we 
recover the ideal case of an ordinary metal with rs,p = ±1 
and Tgp = 0. The integrand at large distances is a nega- 
tive quantity and so Ec{d) approaches zero from negative 
values. From the previous discussion at small distances 
Ec{d) — > -l-oo, therefore there must be a minimum at 
an intermediate distance < dm < 00 since the func- 
tion must cross the x axis. In the unrealistic case where 
e(0) = I one can check that \rsp\ » \rs\,\rp\ making 
Ec{d) always positive for all distances. In this case there 
is no minimum and the force is always repulsive. By this 
analytical analysis and when 9i = —62 we expect that 
as we increase e(0) from one, a minimum develops at an 
intermediate distance dm- This distance shifts to lower 
values as we increase e(0) until, at e(0) = 00 we recover 
the metallic case where complete attraction occurs. In 
the case where 9i = 62 the signs inside (O change and 
make the logarithm to be negative recovering attraction 
at all intermediate distances (for more details see [j^). 
The consistency of these analytical expectations is con- 
firmed numerically with the results shown in figures [2] 
and [3] 

From these we infer that in order to enhance as much 
as possible the minimum, it is necessary to search for 
a situation where 9 = 61 = —O2- The CED for differ- 
ent values of 9 satisfying this condition are depicted in 
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FIG. 3; Effect of parameter e(0) with fixed 61 = —62 = tt. 
The efltect of increasing e(0) is to develop a minimum, which 
shifts to smaller dm as e(0) is increased. 



coating is to gap the surface states. We have estimated 
the parasitic magnetic forces between the magnetic lay- 
ers following [iJl ■ The dipole-dipole interaction is of the 
order of attoN at distances of 50 nm and the magnetic 
Casimir force 2J| is ~ 1/iV, much smaller than the force 
described here which is of the order of 5pN. The proposed 
effect could also be explored in other magnetodielectric 
materials such as Cr^Os which can be described by a 
higher axion coupling |25j. However, these materials in- 
duce more general magnetoelectric couplings [26| which 
we will consider in a future work. 

We acknowledge F. de Juan, M.A.H. Vozmediano, F. 
Guinea, J. Sabio, and B. Valenzuela for very useful dis- 
cussions and suggestions. A.G.G acknowledges MICINN 
(FIS2008-00124) and A.C. acknowledges EPRSC Science 
and Innovation Award(EP/G035954) for funding. 



figure Fig. [2] (b). Under these circumstances the mini- 
mum is more prominent when 9 = n, i.e. when 9 takes 
its smallest possible value. The general analytical anal- 
ysis, supported by the numerical evaluation for different 
parameters, suggest that a simple experimental set-up 
(Fig. [1]) could switch from complete Casimir attraction 
to a stable quantum levitation regime by reversing the 
magnetization of one of the layers covering one of the 
TI. In this process the system will turn from a symmet- 
ric situation where 9i = 92 resulting in attraction to a 
non-symmetric situation where the optimum condition is 
satisfied {9i = — ^2) and a stable minimum appears. 
To conclude, for this appealing situation to be experi- 
mentally accessible one has to search for realizations of 
uiR where the typical distances between the plates are 
at least of order O.l/xm - 1/^m. For a frequency range 
where the axion lagrangian is valid [l7| . the minimum 
is expected to appear at a position where diagonal and 
off-diagonal terms are similar in magnitude, i.e. length 
scales corresponding to 1 — e(w) ~ ■fo^^'- Low values of 
e(0) (typically less than 10) favor this situation since the 
minimum is realized at larger distances. While the elec- 
tromagnetic parameters of TI are still not well charac- 
terized, a low e(0) could be achieved by using thin films 
or by air injection which will lower the bulk dielectric 
response. For TlBiSe2 we estimate from numerical in- 
tegration that the minimum of the CED appears at a 
distance oi d — O.l^m and with a CED of the same order 
than for the usual metal- vacuum- metal system at 1/xm, 
hence being still experimentally accessible. We must note 
however that this estimation requires a high TMEP value 
{9 ~ IOtt) in order to shift the minimum to observable 
distances. Therefore the proposed effect is on the verge 
of experimental accessibility and should encourage ex- 
perimental efforts to attain full optical characterization 
of TI. We stress here that the only effect of the magnetic 
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Auxiliary Material: 
Tunable Casimir repulsion with three dimensional topological insulators 

Attraction versus repulsion from Fresnel coefficients: 



l.Fresnel coejficients for topological insulators: In order to compute the Casimir energy Ec{d), the Fresnel coeffi- 
cients which define the refiection matrices of a topological insulator have to be calculated. These coefficients relate 
the amplitude of the incident and reflected electric fields for different polarizations. Since the normal components of 
D and B and tangential components of E and H must be continuous along the interface one can obtain the following 
reflection amplitudes by defining H and D as in the main text [H Q : 



(r) 
s 

(r) 



1 

A 



"-1«2X- 



2sgn(02-ei)|a|ni 



2sgn(6'2 - 9i)\a\ni 



7iirt2X- 



(6) 



where a = a{62 — di)/Tr, a is the fine structure constant (a 
rij + ^2 + + f^l^^2X+ and: 



Hi is the refractive index of material i, A 



X± 



(7) 



When a — and after a little algebra, these coefficients reduce to the ordinary Fresnel coefficients for a dielectric- 
dielectric interface Since we are considering a topological insulator- vacuum interface, ni = 1, 0i = 0, 02 = ^ and 
n2(w) = therefore a — aO/iT. By its definition, A is always positive and the off diagonal terms have a an overall 
sign governed by the sign of 9. 

To calculate the Casimir energy stored between the plates one computes: 



Ecid) 



Ah JO 

equation (1) in the main paper, where R; is the 2x2 matrix defined in ©. 



(8) 



2. Analytical existence of a minimum and the role of different parameters: To prove the existence of the minimum it 
is enough to assume that the reflection coefficients are given by the expression ([6]) added to the analytical properties of 
the dielectric function, i.e. finite dielectric permittivity at zero frequency (e(0) < oo) and high frequency transparency: 
e(w) — 1 when uj — ?> oo. For an insulator one can assume the dielectric function: 



em 



1 



(9) 



although the derivation does not depend on the analytical form of the dielectric function as long as it fulfils the 
mentioned conditions. In particular, it holds for more than one oscillator. In addition, and for analytical traceabil- 
ity we might further assume the particular situation where 9i = —92 (where labels 1 and 2 identify the Casimir plates). 

It is straightforward to check that if in the expression for Ec{d) we rescale ^ and ky to contain d this gives an 
overall factor 1 /d'^ and forces the reflection matrices to be evaluated at the rescaled frequency and momenta /d and 
k||/c?. Thus all of the d dependence can be transferred to the reflection matrices and the overall 1/d^ prefactor. The 
expression for the energy turns to be: 



1 - Ri 



d' d 



R2 



d' d 



-Ik, 



(10) 
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where k^^^^ = + fc^ is defined through the rescaled variables (which include c the speed of light in vacuum and lur 
defined in the main text) and Eq = Ahc/ {2'k)'^{wrJ cY where A is the area of the plates and the reflection matrices 
are evaluated at these rescaled variables with a rescaled dielectric function: 



e{iS,/d) = 1 + 



1 



IR C 

ujR d 



(11) 



oo. 



The behaviour of e (i^ ) ensures that the reflection matrices are not singular when evaluated at ^/d ^> and S,/d 
Hence we can conclude that Ec{d — )• 0) ^ ±oo and Ec{d — >■ oo) — 0. 

The way the integrand approaches these limits determines the sign of Ec{d). For instance, if the integrand is positive 
at small distances and negative at large distances, necessarily a minimum exists at an intermediate distance dc- lu 
what follows it will be shown that this is exactly what happens when di — —02 (unless e(0) — 1, where both limits 
are positive and hence long range repulsion is obtainned). The first step is to evaluate the integrand in ([SJ. For the 
situation in which 9i — —62 — the reflection matrices describing both topological insulators are: 



R± 



rs(i^, k||) ±rsp(i^, k 
±rsp(i^,k||) rp(i^,k||) 



Introducing this inside ([5]) we obtain the following integrand: 



(12) 



/ — log det 



R+ R e-2fc 



(r) 



log 



-2k 



(r) 



(2. 



2 2\ 



-4k 



(r) 



(r2 



Notice that the last term although always positive, will play no role in what follows since it is always suppressed 



over the first term (note that 



St ' pj ' sp 



<1). 



In the limit of small distances {d — s> 0) and using the high frequency transparency of the dielectric function it 



easy to show that 



<< r 



sp 



(i = s,p): Notice first that the denominator A defined in ([6]) is common to all terms 



so it cannot play a role on the relative magnitude of the coefficients. We need however to study the behaviour of x- 
at small distances, given by: 



X- 



d' d 



(13) 



Remembering that at small distances (large frequencies) we have transparency, n^i^/d) — e{S,/d) — > 1 then it is 
straightforward to see that X- ~^ hence: 



and 



-a 



2 + + X^ 



2a 



2 + 52 + x+ ■ 



since the first are of order 0{a^) and the second are of order 0{a) (remember that a is proportional to the fine 



structure constant a) we have that 



<< \rs 



{i = s,p). Hence the integrand is positive (since the integrand has 



the form / = ln(l + A) where A> 0) and so Ec{d ^ 0) — > +00. 
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Now we consider the limit of large distances {d — )■ oo). In this limit, the reflection coefficients take the form 



l-e(O) 



1 + e(0) + a2 + y<0)x+ ' 
for the diagonal part (with a similar expression for rp) and 



2a 



' sp 



1 + e(0) + a2 + Ve(0)x+ ' 



for the off diagonal. We have defined the quantity e(0) = 1 + ■ 

In this case, depending on the values of e(0) different behaviours emerge. Notice, that from its definition, e(0) > 1 
and so we can distinguish to extreme limits, one where e(0) = 1 and the other with e(0) >> 1. 



When e(0) — 1 one can see that we return to the previous case since the quantity x- in this limit also goes 
to zero. Hence \rsp\ >> \ri\ (i — s,p) is satisfied for all distances and Ec{d) is always positive .Therefore, using 
the fact that that Ec{d) when d — )> oo and that Ec{d) +oo when d — )> we deduce that there is no mini- 
mum and that the force is always repulsive. This is confirmed by the numerical calculations presented in the main text. 

When e(0) >> 1, we see that the opposite condition, \rg p\ » jr^pl, is satisfied. Even in the worst case when X- 
is smallest, e(^) in is always larger than 2a in Vgp. The integrand at large distances is a negative quantity (since 
the integrand now has the form / = ln(l — B) with B > {)) and so Ec{d) approaches to zero from negative values 
{Ec{d) — > — oo for d — > oo). Since at small distances Ec{d) — > +00 we conclude that there must be a minimum at an 
intermediate distance < dm < 00, since the function must cross the d axis. 

Notice that when e(0) is strictly infinity we recover the case of an ideal metal where rs^p = =Fl £^nd Vsp — 0, with 
attraction at all distances. 



To sum up, by this analytical analysis we expect the following situation to occur (when 61 = —02)- As we increase 
e(0) from one, a minimum develops at an intermediate distance dm- This distance shifts to lower values as we increase 
e(0) until, at e = 00 we recover the ideal metal case where complete attraction occurs. 
In the case where 61 ~ 62 the signs inside / change to give: 



I = log 



1 



-2k), 



(2, 



sp 



The predominant term inside the logarithm is always negative and hence the integrand is always negative. Therefore 
at small distances Ec{d) — > —00 and at large distances Ec{d) — > approaching this limit from negative values, 
recovering attraction at all intermediate distances. 

Finally, it can be checked by analogous methods, that the case where one 9 is zero and the other one is not (dielectric 
- topological insulator case) results in Casimir attraction for all distances. 



Comparison with chiral metamaterials: 

It is instructive to compare our results with the results obtainned by Zhao et al. in [4]. In their work, the relevant 
parameter to obtain repulsion is the chirality which mixes transverse electric and transverse magnetic polarizations 
giving in turn the off-diagonal r^p components. When calculating the Casimir force, by differentiation of the Lifshitz 
formula (Eq. (1) in the main text) they obtain the following integrand (in our notation): 



{tI + rl- 2r%)e-^^'^ - {r% + rgVpfe- 



j _ V s > ■ p sp)^ y sp ' ■ s' pj , 

1 1 „2 o„2 \„-2kd\(„2 i -,■ r- \2„-Akd' ^ ' 



Note that when r^p becomes dominant, J turns negative as discussed in [4], enabling thus the possibility of repulsion 
when Tsp is sufficiently large. In chiral metamaterials this occurs for large chirality at high frequencies (short distances) 
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0.02 (1) d/C 0.04 



FIG. 4: Effect of tfie addition of a high frequency osciUator to the dielectric permittivity (4). The units are the same as in 
the figures of the main text. Green squares represent the high frequency oscillator (eq (11611 with cja ~ 0), blue circles the low 
frequency oscillator (eq (|16p with ujt — 0) and the dashed red line represents the model in (|16p 



hence giving repulsion at short distances. The comparison with the topological insulator case is straightforward. 
Analytical derivation shows that the integrand for opposing signs of the topological magnetoelectric response {9i = 
—62) has the exact same form as the integrand given by (|14p (with properly modified reflection coefficients). In the 
last section we showed that, for high frequencies (short distances) r^p becomes dominant and hence, the topological 
part becomes dominant analogous to the case of large chirality in chiral metamaterials. When both topological 
magnetoelectric terms have equal signs {9i =02) one obtains the same integrand with the important difference of the 
signs in front of the off-diagonal terms: 



J = 



2kd 



-f(r 



-ikd ■ 



(15) 



As shown in the last section, even with dominant r^-p repulsion is not possible at any frequency. Therefore the case 
with opposing 61^2 signs is analogous to a high chirality metamaterial. 

Note as well, that the chirality function is frequency dependent while the topological magnetoelectric response is not, 
since it is a topological contribution and it does not get renormalized. There is always an intrinsic high energy cut-off 
of the order of the lattice spacing but since we are interested in the /zm scale these effects are neglected. Hence, with 
opposing signs for 61^2 there is always a region of repulsion for low enough distances. 

Related to the previously discussed works there are other studies where repulsion is realised in the context of 
dielectric and magnetic anisotropy metamaterials (leaving chirality aside). Anisotropy can mix polarizations 

and thus, under the right circumstances induce repulsion as shown by [6|-|8|. The problem of finding anisotropic 
corrections to the Casimir force in topological insulators is not at all trivial. In the present work, we have focused on 
the effect of the topological magneotelectric response on the Casimir force and so we have not included the effect of 
anisotropy which we leave for future work. 

For a full comparison of the cited works with the results here presented it would be necessary to include anisotropy 
in our study in order to draw clear conclusions. Our work is thus closer to that of chiral metamaterials, and so we 
leave comparison with anisotropic metamaterials for a detailed study of anisotropy. Nevertheless one could expect 
that since anisotropy is responsible in metamaterials for repulsion at certain distances (for some specific range of 
parameters), anisotropy can compete with the axionic response to develop richer behaviours. 
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Effect of including two oscillators: 

In the main text we have modelled TlBiSe2 with a single oscillator given experimentally by However in general 
other oscillators at different frequencies can be present when considering other topological insulators. Nevertheless, as 
stated above, the proof for the existence of the minimum still holds even when more oscillators are considered in (4) , 
as long as high frequency transparency and finite zero frequency dielectric response is respected, which is in general 
true for insulators at low temperatures. The presence of other resonance frequencies can modify the position of the 
minimum just as discussed in the main text. To explore this issue a little further we have studied the case when the 
dielectric function is given by: 



in units of cji the frequency of the first oscillator. The parameter /3 = ^ measures the position of the resonance of 
the second oscillator with respect to the first. To illustrate the effect of adding a high frequency oscillator we study the 
particular case where /3 ~ 10^ and = = 0.5 shown in Figure H) We have plotted the cases where both oscillators 
are included red dashed curve) together with the two isolated oscillators (green squares for Wa = w;, ^ and blue 
circles for Wa 7^ = 0). Figure S] shows that the effect of a high frequency oscillator is suppressed whenever Wa 7^ 
and it is only when LOa — Q that it plays a role shifting the repulsive zone to higher distances, as expected from the 
analysis in the main text. Although the minimum still exists in this last case, the dielectric function is so close to one 
that nearly complete repulsion is obtainned, as discussed in the first section. These analysis can be summarized up 
as follows: 1) The minimum still exists when more oscillators are considered, in agreement with the analytical proof 
and 2) The lower frequency will dominate the position of the minimum since its contribution to the dielectric function 
wiU be less suppressed. 
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